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ABSTRACT 


The behavior of a fully-submerged hydrofoil has been 
studied by using Root Locus technique. A multivariable con- 
cept was used in designing the compensators of the system. 

A digital computer simulation program using Digital Simula- 
tion Language or DSL is produced to study the performance on 
the digital computer IBM system 360/67. The Highpoint 


PC(H)-1 is used as the model of this simulation. 
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I. INTRODUCTION 


A. FOIL TYPES 

The main idea of the hydrofoil craft is to reduce the 
resistance of the craft by reducing the area in contact with 
the water and thus increase the speed for the same power. 
One way to accomplish this is fitting foils below the hull 
of a craft, so that the craft is lifted clear of the water 
at high speed. There are two principal foil arrangements in 
use today, classed as surface piercing foil and fully-sub- 


merged foil, as shown in Figure 1. 


£ 


E DX 


ourface-Piercing Foil Fully-Submerged Foil 





FIGURE 1 


TYPEOF"AYDROFOTL CRAFT 


The surface-piercing foil is supported by vertical 
struts and the lifting surface extends in a V-shape from the 
water surface down through the water to the point of conver- 


gence of the V. The surface-piercing foil provides a 





self-stabilization for the craft in pitch, height and roll. 
The self-stabilizing aspect of the surface piercing foil is 
easy to see. Suppose the craft pitches up slightly, then 
less area of the forward foil and more area of the after foil 
is submerged. This results in less lift from the forward 
foil and more lift from the after foil so the craft will 
settle back again, the same stabilization occurs in roll. 
For the height as the craft speed increases the whole craft 
rises up until the amount of wetted surface of the foils is 
just sufficient to create a lift equal to the weight of the 
crar 

The fully-submerged foil cannot adjust itself to chang- 
ing lift demands by a change in area, but only by a change 
in incidence, either of the whole foil or of attached flaps. 
This involves the use of some sensing arrangement to control 
the incidence by monitoring the height of the craft above 
water surface. The fully-submerged foil is a more sophisti- 
cated design, and the control mechanism is likely to be com- 
plicated. But the fully-submerged foil has a very distinct 
advantage when operating in heavy seas. The surface-piercing 
foil will respond to even relatively small changes in the 
surface of the sea and give a very rough ride. The fully- 
submerged foil, on the other hand, is much less affected by 
the change in the surface of the sea. This paper will deal 


only with the fully-submerged foil craft. 





< 


Le 





B. FOIL ARRANGEMENT 

Most of the fully-submerged foil craft have two foils, 
one at the forward and the other at the aft. Both foils are 
supported by the vertical struts. Figure 2 shows the con- 
figuration of main foil and secondary foil; Figures 3A and 
3B show foils in canard arrangement and in airplane arrange- 
ment respectively. 

The simulation in this paper is based on the canard ar- 
rangement, which has a secondary foil as a forward foil and 
a main foil as the aft foil. The HIGH POINT PC(H)-1 was the 
model for this simulation. Drawing of the craft and the 


foils are shown in Figures 4A and 4B. 
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FIGURE 2 


FOILS CONFIGURATION 
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FIGURE 3A 


FOIL IN CANARD ARRANGEMENT 





FIGURE 3B 


FOIL IN AIRPLANE ARRANGEMENT 


18 





(quvoguvis auv 
-LDHUOd)NONSTIV 


AIVI YNS TOULNOO* ha y 


NOI IVENDIINOO INUIS ANV TIQ 40 ONI AVEG 


Vig 290901 


(1VOIWISWASY 
Ad HIE 'IOHLNOO QUYMIOJ 


(NMOHS LON SI ‘some 
ATLOuHIq'usaanu zavds) 
udqans Soda ONITIVAL 


1O 








WAS CAYMIOT 


NOILVYUMDIANOO AVI GNV TIOI JO PNIMVYC 


ay NII 


IF 
€ 


, 


E ea 


CUM  ENIATA — 


` 


^v 


or 


A 


INC SYB 


7107 LIV 








20 





II. HYDROFOIL EQUATIONS STUDY 


A. COORDINATE SYSTEMS 

Three coordinate systems will be used in this study of 
hydrofoil craft dynamics. They are: the Body Axes, the 
Water Axes and the Earth Axes. Each is in a right hand 
orthogonal system. | 

The body axes coordinate system (X,Y,2) shown in Figure 
S has its origin at the center of gravity of the craft and 
is fixed relative to the craft. This system is initially 
coincident with the origin of the earth axes system, however, 
it moves with the craft as time progresses. X is positive 
forward, Y is positive to starboard and Z is positive down- 


t 


ward. 






FOSITIVE p 
(ROLLING ro stED) 


X - AXIS 


POSITIVE Y 


POSITIVE R 





zo STED,) 


FIGURE 5 
ORIENTATION OF THE BODY AXES WITH RESPECT TO THE 


CRAFT AND DIRECTIONS OF POSITIVE VELOCITILS 
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The earth axes coordinate system (XL, Y Zr) is fixed 


E? 


relative to the earth surface. The origin may be chosen at 
any convenient location, provided the body system origin is 


initially at the same point. The Xp axis lies in the hori- 


zontal plane; it 1S initially coincident with, and positive 


in the same direction as the body X axis. The Yt axis lies 


in the horizontal plane and is positive to starboard when 


E direction. The Zn 


axis is normal to the horizontal and is positive downward. 


the observer is facing the positive X 


Water surface motion and craft motion relative to calm water 
are described in this coordinate system. 

The water axes coordinate system is aligned with the 
relative velocity vector and resolves into "Lift," "Drag," 


and Side-Force" direction. 


LIET 


E 
SIDE FORCE 


D 
DRAG 


It will be specified that the order of rotation from 
water axes to body axes is: 

1) The (D,S) axes rotate about the L-axis through the 
angle ß to the position (D,,5,) where Si is coincident with 
the body Y-axis. 


2) The (D, ,L) axes rotate about the S,-axis, or more 


I 


pertinently about the Y-axis, through the angle a to the 


position (Do, L4) which corresponds to the (-X,-Z) axes. 





B. EQUATION OF MOTION 

The development of the six equations of motion for a 
hydrofoil follows the same procedure as those for a displace- 
ment type hull. This development is clearly outlined in 
Ref. 2 and 8 and therefore will not be done here. 


The equations of motion are: 


Translation 
O ou, (2-1) 
m X 
en (2-2) 
m Y 
O E (2-3) 
m A 
Rotation 
I «Be l : 
P = I [L-QR( I, - I, ) t CR*QP)I,,] (2-4) 
XX 
|a Wl 7 _ Bee _ 
Q = = RP(Iyy-1,,)-(P2-R2)I,] (2-5) 
ST “Ta " " ° 
R = O Ty) (QR-P)I..., ] (2-6) 


The terms RV, PV, and PW are very small and may be neg- 


lected. Similarly, the product of inertia quantity, lez is 


less than 10% of I and less than 4% of I Consequently, 


XA 22° 


when considered in conjunction with the magnitudes of (R+PQ), 


e R and (QR-P), the terms containing I,, may be neglected 


XZ 
Tor this craft. 
Equations (2-1) through (2-6) are linear and rotatuonal 


accelerations in body axes. The velocity terms in body axes 
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can be obtained directly by integrating each equation one 


time. 
U = / Udt (2-7) 
e vai (2-8) 
W = f Wdt (2-9) 
p = f Pat (2-10) 
ONSE Odt (2-11) 
R - / Rdt (221095 


Linear displacements and angular positions with respect 
to body axes could be obtained by an additional integration, 
but such quantities are of limited usefulness. The dis- 
placements and angle usually of interest are those with 
respect to earth axes. In order to calculate the craft's 
position in earth axes coordinates, the craft velocity vec- 
tor must be transformed from body axes coordinates into 


carthn#l-ies#coordinates.2s follows 


Ug = Ucos@cosy+V[coswsin®siny-sinycosg] 


+W[cosbsindcosd+tsinysind] (2-13) 


VE = Ucos@siny+V[cosywcosd+sinysingsinp] 


+W[sinywsin@cosg-cosysind ] (2-14) 


Wn = -Usind+Vcosdsind+tWcosdcosd (2-15) 


Only the craft's velocity in the downward, Lo direction 
is@required tofcalcullate its height above the water surface. 


That is, the vertical position in earth axes follows from 


Ze dt (2-16) 
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The craft's position coordinates are obtained from 


Ze f Updt (2-17) 


Ze JS V„at (2-18) 


The Euler angles, 0, 9, and Y are also of interest and 
can be found by integration of Euler rates Ó, o. and y. The 
rates in turn must be derived from the body angular rates. 
The Euler rates are not easily obtained because they occur 
about axes which are not orthogonal. This fact can be ap- 
preciated by referring to Figure 6 and recalling how the 
angles were defined. w is a rotation about the Zo axis, 9 
is a rotation about the Y, axis, and 0 is a rotation about 


the X axis. They result in the following set of equations: 


4 = P+ysin? (OD 
ô = Qcosd-Rsind (2-10) 
bh = (Qsinòg+Rcosò)cos9+(d-P)sin9 (OE 


Studies of this craft show that even under the most 
drastic abnormal situations, pitch angles are not expected 
to reach 10% and under normal operating conditions, roll 
angles will not exceed 3°. Therefore, the following assump- 


tions are justified; 


sind O 


(oor 


lI 


così T 


The Euler angles are obtained by integration of equations 


(2-19) through (2-21). 
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FIGURE 6 


DEFINITION OF POSITIVE DIRECTIONS OF EULER ANGLE 





C. CALCULATION OF VELOCITY COMPONENTS 


The total velocity components of each foil i 


n body axes 


can be expressed in terms of craft linear and angular veloc- 


ities by considering the craft geometry as re 


Figure 7 and Table LE 


The resultant equations for the example s 
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presented in 


(2-23) 
(2-24) 


(2-25) 


(2-26) 
(2-27) 


(2-28) 


(2-29) 
(2-30) 


(231 


(2222) 
(2225) 


(2-34) 











FIGURE 7 


DEFINITION OF SYMBOLOGY 
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and L are variables and in general will not be 


ZOS 


equal. 
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a | 


These foil velocities are in body axis terms and can be 
used directly to calculate the angles of attack and angles 


of side slip. 


D. ANGLE OF ATTACK AND ANGLE OF SIDE SLIP 

The angle of attack of the foil is shown in Figure 8A 
and defined as the angle whose tangent is the total relative 
velocity between the water and the foil in the body Z direc- 
tion, divided by the total relative velocity between the 


water and foil in the body X direction 


W 


o = arc tan = (2-35) 
E 


If the foil is not aligned with the body X-axis, then 
Pie com lies a isecd angle of attack and the total angle of 


attack becomes 


“TOTAL 7 ČFIXED T| Š (2-36) 


The side slip angle is shown in Figure 8B and can be 


represented by 


V. 
CI 
p Sare Van ee (2-37) 
i /U2 ..+V2 +W2. 
viva P pool 
During actual craft operations, the angles will be less 
than 10°, and will permit use of the small angle approxima- 
tion. Also the V and W velocities can be considered very 


small in relation to the U velocity. Applying these simpli- 


fications, the equations for a and 6 become: 
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FIGURE 8B 


ANGLE CF SIDE SLIP 
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E. EXPANSION OF FORCE AND MOMENT TERMS 

With the aid of Figure 7 and Table I, the forces and 
moments are set down. It is important that when a numerical 
value is inserted for one of the alphabetical dimensions, 
the sign of the value must be positive or negative according 
to positive and negative directions from the axis system 
origin at the C.G. Also note that horizontal foils make no 
Conia puton to thefbody Y force and vertical struts do not 
contribute to the body Z force. One of the assumptions made 
upon commencing simulation was thai the struts contributed 
LOL LL EAN hat the foils contributed zero side force. 
Therefore, all body Y motion originates at the struts and 


all body Z motion originates at the foils. 


The force and moment equations become: 
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All of these forces have a hydrodynamic origin except 
for the thrust and gravity terms. The thrust is always 
associated with the body axis and requires no further expan- 
sion. The gravity terms can be simplified using small angle 
approximation for the Euler angles (see equations A-31,32,33). 

The general expression for hydrodynamic forces in water 


axes coordinates is 


F, = qAC, i = L,D,S (2-44) 


where L,D,S represent lift, drag and side force respectively. 
The lift, drag, and side force coefficients very as func- 
tions of angle of attack, angle of side slip, submergence, 
velocity, flap elevator and rudder positions. 

Ideally, mathematical expressions would have been de- 
veloped to correctly depict the interrelation of all the 
variables which affect the force coefficients. However, 
hydrofoil technology has not advanced to a position which 
would yield such expressions. Some of the major problems 
encountered in deriving a mathematical expression are: 

1) The occurrence of cavitation and ventilation which 
are completely unpredictable. 

2) Even more unpredictable is the cessation of cavita- 
tion and ventilation. 

3) The nonlinearity of the hydrodynamic coefficients. 

4) Lack of sufficient test data to completely describe 


all of the above. 








F. SPECIFICATION OF UNITS 

Considerable confusion resulted from the interpretation 
of the various references concerning the units used in de- 
fining Euler angles and ER rotations. This problem has 
been alleviated here by the specific definitions in the 


Table of Symbols. 





A. SIMULATION OBJECTIVES 


HIT. 


HYDROFOIL SIMULATION 


The object of the simulation was to commence the runs 


with the craft foilborne and in a steady state condition. 


Once this condition was achieved, the step response of the 


craft would be studied by applying step functions to excite 


motion along and about each body axis. That is, perturba- 


tions should be introduced to separately excite: 


1) 


2) 


3). 


4) 


9) 


6) 


Motion 
Motion 
Motion 
Motion 
Motion 


Motion 


along 
about 
along 
about 
along 


about 


the 
the 
the 
the 
the 


the 


X-axis 


x-axis 


Y-axis 


Y-axis 


Z-axis 


Z-axis 


After the step response was obtained, a sinusoidal sea 


would be inserted to observe the craft's motion in a seaway. 


The following assumptions were made prior to beginning simu- 


lation: 


1) 


2) 


3) 


4) 


Craft equations of motion are valid only for the 


foilborne mode. 


Weight of tbe craft remains constnat. 


Hydrodynamic coefficients are based on fully wetted 


Suri aces LE 


mo Cavitatlion or ventilation. 


NOorconstrarint was placed on craft's heading. 





5) Foils are considered flat surfaces vice dihedral/ 
anhedral. 


6) The craft may be operated in calm water or seaway. 


B. COMPUTER SIMULATION 

The computer language used in this simulation was Digital 
Simulation Language or DSL, and based on the computer pro- 
gram in Ref. 1. 

Before the main program could be assembled, subroutines 
had to be written to find the hydrodynamic coefficients of 
foil lift and drag, and strut drag and side slip. Using the 
Simulation curves in Ref. 3 as a guide, the curves were 
stored as data points in the setup routine. The subroutines 
"INTERP" and "INTERPL" are interpolation routines to obtain 
the proper coefficient from the curves. For INTERP, values 
of angle of attack and flap angle are used to enter the 
curves and obtain a value of foil drag. Angle of side slip 
and submergence are used to enter INTERP1 to obtain strut 
drag. The curve for the lift coefficient of a foil in the 
fully wetted region is straight line so this coefficient is 
una by merely solving the equation for the straizght line. 
Curves plotted from the subroutines are shown in Appendix B. 
Listings of the subroutines are shown in the computer pro- 
gram section. The main program was then assembled using 
Figure 9 as a reference for data flow. One part before the 
end of the main program was opened for an automatic control 
block. A complete listing of the main program is shown in 


the computer program section. 
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C. AUTOMATIC CONTROL 

Control signals originate at the pilot's controls, at the 
motion sensors, and at the position transducer. These are | 
combined and processed in the electronics to produce the 
signal for the servo-valve. The pilot will be able to move 
the helm, a lever, or a knob to introduce steering, altitude 
and attitude commands, the motion sensors sense errors be- 
tween the commanded and actual craft motions and produce 
electrical signals proportional to these errors. The posi- 
tion transducer provides the feedback signal so that the 
control surface stope moving when it has reached the position 
corresponding to the processed error signals. 

It is necessary to consider how the control system com- 
ponents aiiect tae equations of mctior through the controi 
surface deflections. In general, a signal originating at 
the pilot or at the motion sensors will be modified by the 
individual dynamics of the components in the signal path be- 
fore the control surface motion actually occurs. Motion 
sensors to be discussed are vertical gyros, rate gyros, ac- 
celeromaters and height sensors. 

me Vertical" Gyros 

NS? io EE EE eet reste will 
produce one signal directly proportional to the pitch angle 
Oi the cralt und one signal directly proportional to the 


rose allele Ge the crait. 





2. Rate Gyros 


A rate gyro will sense an angular rate about one of 
the three body axes. (Three-gyro "rate packages" are avail- 
able if required for a particular application.) The dynam- 
ics of the gyro are such as to produce an output signal 
which is related to the actual angular rate by a second order 
differential equation. Using the roll rate P, as an example, 


the equation is 


ke, d 


ap 2 
+ — _ 
250. it 4 wP) (3-1) 





2 


K E = wi 
5 7 E 7 S u 54220 Stw (3-2) 
where | 
P,Q,R = Outputs of rate gyros (RAD/SEC) 
P,Q,R = Craft rotational rates (RAD/SEC) 
a = Natural frequency of a particular gyro 
(RAD/SEC) 
E = Damping ratio of a particular gyro 
(Dimensionless) 
S = Laplace transform variable 


3. Accelerometer 
Outputs of accelerometers are also assumed to relate 
to the craft motions by second order differential equations. 
In general form, the equation is: 


da, dä, 
a. - Li, 204 —i t waan) (3-3) 
1 "e dt? n dt Ce? 
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where i will become a symbol corresponding to the location 
of a particular accelerometer. 

To identify the accelerations sensed by the instru- 
ments, it will be assumed that there is one accelerometer at 
the C.G. mounted such as to measure lateral acceleration, 
and one accelerometer above each foil mounted in the X-Y 
plane so as to measure accelerations parallel to the Z-axis. 

The lateral acceleration at the C.G. is: 

F 


= Y E 
a, = (3-5) 


Taking angular acceleration into account, the accei- 
erations in the Z-direction at points above a center, port, 


and starboard foil are: 


La 


“eo + Lon" (3-6) 
A St ep aw 
Pen XPF YPF 
E ° : | 
E Rr Sue Sr «oor 


4. Height Sensor 
Ultrasonic and sonar types are either in use or ex- 
pected to be used. For this study it will be assumed that 
the dynamic characteristics of the sensor, whatever the 
type, will be such as to give a direct indication of the in- 


stantaneous height of the sensor above the water. 
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IV. TRANSFER FUNCTION AND ROOT LOCUS DESIGN 


A. PITCH-HEAVE MOTION TRANSFER FUNCTION 

The best way to study the stability of the craft was to 
find the transfer functions of the craft and use either Root 
Locus plots or Bode plots to be an aid in a stability study 
and also in designing the compensators of the system. 

The transfer functions of the craft can be derived di- 
rectly from the equations of motion. In the pitch-heave 
mode, equations (2-3) and (2-5) were the equations that de- 
scribed the motion in the Z-axis and motion about the Y-axis. 
After neglecting the small terms, equations (2-3) and (2-5) 


were as follow 


- 


wet, : 

W = = FZ'Q U (3-9) 

c.i 

Q = Te M (3-10) 
YY 


From equations (2-35) and (2-36) the angle of attack of 


each foil can be written as follows: 


ag = FS + ALL (3-11) 
a g = Rs + Alig (3-12) 
Gp = “RP + AL, (3-13) 
“y “RI + AL (3-14) 
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Wan = W + Lupp P - Lepp Q (3-15) 


Was = Wt Lysr'P - Lygp € (3-16) 
Wen © W + Lxgp 2 (3-17) 
Wap = W= Loop" 0 | (3-18) 


Assume there is no significant coupling from the roll, 
yaw, sway or surge motions. Equations (3-15), (3-16), (3-17) 


and (3-18) will be 


Mr sss pic 
MES sr 9 
ibe ees ~ 
"cM NE cri: 


And equations (3-11), (3-12), (3-13) and (3-14) will be 
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(3-19) 


(3-20) 


(3-21) 


(3-22) 


Let ALL A,» ALS = AL, = AL, — 
“XCE. g Qi 
XSF. Q = Q2 
oc = W -Qi*Ai 
Ag = W -Q,4A, - dp = Du 

From Equation (2-40) E, = cr zp srt our! ez 

Fzcr ^ "tpoz ?c "rez 
Topp " "Fppz'"p"frpz 
Fzsr "© "Tpsz "s sz | 
Four © © pme “7 LMZ 

Equation (2-44) shows the general expression for hydro- 


dynamic forces, 


hydrodynamic 


Chop 


CLop 


Chan 


Chirp 


and from Ref. 


coefficients are: 


Dye ocn 


4.10 a 


= 4,10 a 


Il 


5.15 a 


C 


p 


S 


M 


+ 0.036 Sp 1007 183 
+ 0.018 Sp 


+ 0.018 Ss 


+ 0.036 SA 


3 the linear approximation of 





Chop = 0.012 + 0.0015 6, 
Cha = 0.02 + 0.0028 8, 
Cop = 0.02 + 0.0028 8, 
Chg = 0.02 + 0.0028 64 
Let 
0.5*p*U* * Ac, =D, 
m. = 0.5*p* Uode = Dr 
EE, em =, 


After dropping the nonlinear terms 


Poop = 7(5-162D,)W +(5.162D,)Q,-D, (0. 0015A, +0. 036)67 


-D, (5.1624,+0.13) 


Fypp 7 -(4.12D, )W *(4.12D, )Q, -4. 12D, A, 

Fogg ^ -(4.12D, )W *(4. 12D, )Q, -4. 12D, A, 

Foye = -(5.17D,)W *(5.17D, )Q,-D, (0. 00284, *0.036)6 ,-5. 17D, A, 

FK,  - -(5.162D,48.24D, *5.17D, ) *(214.69D, -139.92D, -87. 79D, )3 
-D, (0. 0015A, ^0. 036)6 ,-D, (0.00284, *0. 036)6 , 
-[D, (5.1624, 40.13) *8. 24D, A, *5. 17D, A, ] 

5 = _(8.162D,+8.24D,+5.17D,) y „ ((214.69D,-139.92D,-87. 79D) 

mU mU 

Meda -G ; - D1(0.0015,50.000), _ Da (0.00284, +0. 036), 


.. [D1 (5. 162A, *0. 13) *8. 24D, A; t5. 17D4A;] 


= (3-23) 


44 





M 
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Locs'txpsLzpPs 


F L 


FrsrbzsrtFxurbzur'"xcs 


E L L 


-M AFT *xcr*zcF* 
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FWD XSF 


tFesslzsstTy zT 


Pocrlxcr 


F )L 


Fase} ZMF 


Si XMF 


ZPF' 
[(214.69D, )W -(214.69D,)Q,+(0.212D,A,+1.497D,)8, 
+(214-69D,A,+5.41D, )]+[-(139.92D,)W +(139.92D,)Q» 
-139.92D2A2-(87.79D3)W +(87.79D3)Q2-(0.048D3A2+0.611D3)8, 
-87.79D3A421+[(281.1DA2)W -(281.1D2A2)Q2+(140.55D2A? 
-0.686D>)1+[(159.65D,A,+2.02D, )W -(159.65D,A,+2.02D,)Q, 
+(0.558D,A1-0.023D,)8,+(79.83D,A1+2.0241D,-0.186)] 


*E(176.5D3A2)W -(176.5D3A2)Q27(0.62D542-0.048D5 6, 


*(88.27D442-0.343)]-14. 29U^-17. 7T. 


[D, (216. 71*159. 654)-D; (139. 92-281. 1A; )-D4(87. 79-179. 5A, ) ]y. 
1 U | 


YY 


_ {D1 (8928 . 96+6639. 84A, )t+D, (2375. 84-4773. 08A, )+D3 (1490. 67 
I U 
vid 


—2996.97A, )] D,(1.474*0. 7741) D,(0.659-0. 572A, ) 
YY YY A 


+ [D, 5.41{A, (216. 71+79.83A, }}-D:{0.686+A, (139.92 


-140.548A2 )}-D3 A, (87. 79-88. 27A, )-0. 529-14. 29U2 +17. 77] 


yy 


(3-24) 





From (3-23) and (3-24) 


EA AAA EA EE 


E A 


-C.5 ES C0 + 7 + SR = -Ca ò -Ca ô - U, 


F A 


By taking the Laplace Transform of the two equations 


above and rearranging into a matrix form 


ENS t6) -C18 9 C3 -C, || 6 Ú, 


=(C6S+C7) S(StCs) Z -Cg -Ca |} ô -U; 


The characteristic equation 


A = S[S{S°+(C2+Cs}-C1C6S-C1C7] 


where: 
c. = 1D1(216.71 139.65A,)-D,(139.92-281.1A,)-D,(87.79-176.5A,)] 
: Tu 
YY 
c, = [22 (8928. 96+6639.84A, )+D, (2375.84-4773.08A,) 
: ës LU 


+D:(1490.67-2996.97A; )] 


[D1(1.474+0.77A,)] 


Cs = 
Lyy 
[D4(0.659-0. 5724; )] 
Cu : 
YY 
c. - [5-162D, +8.24D,+5.17D] 


mU 


om = [ (214. 69D,-139.92D>-87.79D3) +U] 
mU 
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76 
m 


D,(0.0015A,+0. 036) ] 


Cs = m 


After inserting the value of parameters, a speed of 36 
knots was chosen to be used in this simulation to insure 
that cavitation did not exist, and it was assumed that there 


was no external disturbance. The transfer function matrix 


becomes: 
S(S+10.791) 20.0635 || 0.083 -0.0394 sa | 
1 A E || 7 -1.04  -1.215 ||6, 


and the characteristic equation was 
A = S[S*+18.4S°+77.31S-2.03] 


= S(S-0.0261)(S+6.55)(S+11.876) 









and 
E P _ 0.083(S+6.82) 
€ (S*0.0261)(S*6.55)(S*11.876) 
de S(S+7.61) 
A 
S(S+10.731) 0.083 
ES _ = 1.04(S-0.495)(S+5.19) 
Me S(S-0.0261)(S*6.55)(S*11.876) 
Sa O 76.338432.2) 





ag 





-0.0394 -0. 0635 


0 -0.0394(849.55 
SA (S-0.0261)(S+6.55)(S+11.876) 
6570 | 21.215 S(S+7.61) 
A 


S(S+10.791) 0.0394 


E e _ _1.215(S+0.079)(S+13.186) 
€, 3(S-0.0261)(8+6.55)(S+11.876) 
S070 (1 (76.338+32.2) 1.215 





A 


Ya — -Usin0-4Vcos0sinotWcos0coso 


V, the velocity in Y-direction is equal to zero, and 


angle 6 is assumed small. 


Sin 0 = 0 E Cos 8 = 1. 
= —3 o Á 
WE U-90 + W 
E de 0. 
-heg 7 T/k Ga 


The height at the sensor, 























Sm -heg Lon" 
MEN: 
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Amos MS EI Row 
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The four transfer functions of the craft are: 








0 0.083(S+6. 82) 
Sp (S-0.0261)(S+6.55)(S+11. 876) 
6 ,=0 
-hs e OMG ° (Sr0”2 (C Sr6752) 
F S(S-0.0261)(S+6.55)(S+11.876) 
§ ,=0 
A 
Fo" = -0.0394(S+9.55) 
J (S-0.0261)(S+6.55)(S+11.876) 
5_=0 
F 
. _ -1.464(S+2.72)(S+5. 348) 
di S(S-0.0261)(S+6.55)(S+11.876) 
$_=0 
F 


The characteristic equation shows that this system has 
one pole in the right half plane, and from the Root Locus 
plots three of the plant transfer functions have roots in 
right half plane, only the effect of forward foil to the 
height sensor has complex roots in the left halp plane but 
very close to the imaginary axis, which was only marginally 
stable. From this point, the stability of the craft without 
automatic control system was questioned, and led to the con- 


clusion that the craft must have an automatic control system. 


B. ROOT LOCUS DESIGN OF FORWARD FOIL AUTOMATIC CONTROL 
os clei No Bona themrorward foil as a control surface, 

the aft foil was fixed at the angle of 0.074269 radian [Ref. 

1]. By closing the height loop and the pitch angle loop as 


in Figure 14 and increasing gain K of the system to 5, from 
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Figure 10 the roots of the height loop moved further into 
the left half plane and became repeated real roots at -1.6. 
This value of gain K pushes the root of the pitch angle loop 
(Figure 11) just across the imaginary axis to the left half 
plane. The system becomes stable. The response to the 
change in height command of 2.0 feet was shown in Figure 15. 


The craft went to steady state in 5.5 Sec. 


Léi d; ES 
FE "e mo de use 
N qa = PLANT 
Cher (SA li i a 
2 | 
+ re 
FIGURE 14 


BLOCK DIAGRAM SHOWING FORWARD FOIL AUTOMATIC CONTROL 


C. ROOT LOCUS DESIGN OF AFT FOIL AUTOMATIC CONTROL 
in like manner, if the aft foil is used as a control 
surface and the forward foil is fixed at the angle of 


0.03711 radian, the system is unstable by itself, because 
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STEP RESPONSE OF THE CRAFT 


(FORWARD FOIL ONLY) 
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the open loop roots of both height and pitch angle were in 
the right half plane. From the Root Locus plots Figures 12, 
13, positive gain will have the roots further into the right 
half plane. After closing the loops and using gain K equal 
to -47, the height loop had complex roots at -31j3. The 
right half plane root of the pitch angle loop moved toward 
the left half plane, across the imaginary axis to the point 


-0.2. The system becomes stable. (See Figure 16.) 


Azer X) 
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FIGURE 16 


BLOCK DIAGRAM SHOWING AFT FOIL AUTOMATIC CONTROL 


Figure 17 shows the response of the craft to the change 
in height command of two feet. Because the flap angle was 
dimitea at 25°, and the gain of the aft foil was too high, 
in the first five seconds the flap exceeded its limit and 
the flap angle was set to 5° which had the same result as 


reducing the gain of the system. Then in the first five 
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FIGURE 17 


STEP RESPONSE OF THE CRAFT 


(AFT. FOIL ONLY, INPUT = 20) 





seconds the craft behaved like a very lightly damped system, 
but after five seconds the craft; became more stable and 
reached steady state in eight ES Figure 18 shows the 
response of the craft when the change in height command was 
0.5 foot. In the first 1.0 sec. the flap angle still ex- 
ceeded the limit, but the craft went to steady state in 3.0 
seconds. 

Parts B and C have shown that the stability of the craft 
in pitch-heave mode can be achieved by using either forward 
or aft foil as a control surface with the other foil fixed. 
But to make the craft completely stable in six degrees of 
freedom both forward and aft foil have to be used as control 
surfaces. The dynamics of the system will produce a cross 
coupling effect between the two inputs, the systems now have 
multiple inputs and multiple outputs, which is a classical 


multivariable system and will be studied in the next chapter. 
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FIGURE 18 
STEP RESPONSE OF THE CRAFT 


(AFT. FOIL ONLY, INPUT = 0.5) 
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V. MULTIVARIABLE SYSTEM 


A. BASIC CONCEPT 

A system which has multiple inputs-multiple outputs is 
much more complicated to design than a single input-single 
output system, because a change in one input may have an ef- 
fect not only on the corresponding output but also on the 
other outputs. In designing this kind of system, special 
techniques are needed. Ref. 9 has all of the theories and 
procedures to work on this kind of problem, and this paper 
will follow the procedures in that reference. 

The idea was that, the plant has a nxm transfer function 
matrix and was compensated by a mxn diagonal matrix compen- 
sator with unity feedback. The configuration of the compen- 
sated system is shown in Figure 19. 

Consider the hydrofoil craft in pitch-heave made only, 
and not coupled to the other modes. The craft will have a 


2x2 transfer function matrix. 
EN EC 
G = (4-1) 


g g 
P21 Da: 


If the diagonal compensator matrix 


0 
SE 


SE (4-2) 


C22 


60 





is used, the system configuration becomes 





FIGURE 19 


BLOCK DIAGRAM SHOWING A MULTIVARIABLE SYSTEM 
WITH DIAGONAL MATRIX CASCADE COMPENSATOR 


pii ce 012022 


E > (4-3) 
| “pie, uo A 


1. The Characteristic Equation 
The stability of a multivariable system is determined 


by the zeros of N,(s) and N(s), which are defined as follows: 


a A det oe (4-4) 
x AS) CS) 
N(s) A A (4-5) 


and As), a represent the characteristic polynomials of 
the rational transfer function matrices G. (s) and TS 
respectlvely. 

If 


1) Cancellations are selected systematically, as in 


Ref. 9, 
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2) Poles of G. are carefully selected, 
then (4-4) alone determines the stability of the system. 


Thus 


det(I+G G ) = O 
INTE 


will be referred to as the characteristic equation of multi- 
variable system. 
2. Connection Between Single-Loop and Multivariable 


Cases for 2x2 Plant Ss: and 2x2 Diagonal G 


det(I+G G = 1+ + +(det-G 
p e) pra Ci) IS ec p^8c, ,8c,, 
g +(det-G_)g 
= (1+ e DE Ci 
( 80,4, ,? 1+g g SEN 
Pri Cii 
(4-6) 
I = (1+ 1+G 4-7 
Mo ee, (4-7) 
where AE 
ES O 
g +(det-G Eo Pew B c 
G ni 
e 
E oo 511 C11 


The roots of the characteristic equation are the 


zeros of the rational function (1-4 1+G È 
( Be, Be, Il ende, >? 


From (4-6), can be seen that the roots of 


det-G 
1 + A fn = 0 (4-8) 
g C11 
P22 
and 
ie = 4-9 
Spunea Š ( ) 
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constitute part of the zeros and poles of x respectively 
det*G 


and and g are known functions. 
D22 P11 

Therefore, any single-loop method can be used in de- 
signing EE to place the roots of (4-8) and (4-9) at desir- 
able locations, and also has toform a reasonably good 
pole-zero pattern for EE After GER is designed, all 
poles and zeros of E are known and E can be designed 
from 

M = 0 (4-10) 
B. ROOT LOCUS DESIGN FOR PITCH-HEAVE MODE 


From the craft transfer function matrix 


(hg Eg 
Sp SA 
G = 
p 
e e. | 
6. SA | 
_ Ths _ __6.68(8+0.72)(5+6.52) 
Sp: $. — S(S-0.0261)(8*6.55)(S*11.876) 
_ “5g _ -1.464(5+2.72)(S+5.348) 
Epi. 5,  S(S-0. 0261) (S+6.55)(S+11. 876) 
MOS 0.083(S+6.82) 
pi ôn (8-0.0261)(5+6.55)(5+11.876) 
- UD x -0.0394(849.55) ` 
E 84 — (8-0.0261)(8*6. 55) (S* 11.876) 
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det °G g 


p P11%p22 OP 2221 


-0.141(S-0.02325)(S+6.51)(S+12.02) 
~ §(S-0.0261)* (S+6. 55)*(S+11.876)? 


det*G "e m a 
p 5 57/90590902925) (STO. 51)(5T12.02J9 


ED.) S(S-0.0261)(S+6.55)(S+11.876) 


The design procedure starting from PONE by plotting the 
root locus of equations (4-8) and (4-9) as in Figures 20 and 
21, from Figure 20 can be seen that, as the gain was in- 
creased the root in the right half plane moved toward the 
left half plane, met pole at the origin, broke away and 
circled around the zero at -0.72, with the value of the gain 
equal to 5. Two dominant roots were at -1.25 and -2.2, and 
other two roots at -6.75 and -11.45 and the effect of this 
gain pushes the right half plane root of equation (4-8) to 
the left and was cancelled by the zero at +0.02325 as in 
Figure 21, the pole at the origin moved further into the 
left half plane to the point -2.3 and the other roots were 
in the left half plane at -6.8 + jO.8 and -11.99. Therefore, 
by using the gain equal to 5 as an: the height loop was 
stabilized and, this also gave a good pole-zero pattern for 
Sag: 

Equations (6-30) - (6-35) of Ref. 9 show the procedures 


to find poles and zeros of MEG which were listed below 


G, = it 4-11 
1 CO ( ) 


G2 


+(deteG 4-12 
UE jale pee), ( ) 
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Der = Gm (4-13) 
Pa = Zc *PG - Gg fl Mes i (4-14) 
eq 
Ze = 26, Pa (PG GEN (4-15) 
eq . 
p ER = P +P « (4-16) 
ER 0119011 8011 8011 
Gi l ip c M ( ) 
det * G 
= + — 
46, Zi = u, (4-18) 


And from the system transfer function matrix 


Ke 
Il 
J 
FO 
N 
KJ 
+ 
tul 


- = {0,0.0261,-5.55,-11.876}+P 
G 3 > > 2 
3 EDIT "cin EM c “eni 
Z = Z, + 
G; i SD CM 
DE +P(det-G Ee ee 
G5 ( p^8c 1 g ( p^8c, , ! 


P22 : P22 


0102615 6. 51 S76 10. (0.0261)2,(=6.55)*,(-11.876)*) 


Uu -[{0.0261,-6.55,-11.876}f] ONU D261) 55) 35 
C11 


(-11.876)?)+8, 1 


10560209262, 7 28.5596 ,(-11.876)° T -[{(0.0261), 
C11 


(-6.55),(-11.876)}] 
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Lo, (0.0261)%,(-6.55)%,(-11.876)%)+P, 


C11 
det*G 
x n 
Des =) 801) 
P22 
Pç f! PG, = (0,0.0261,-6.55,-11.876)+P, RPC 
Gia 1 
DU TUE WS +{0,(0.0261)°,(-6.55)°,(-11.876)°}+P 
eq PINE] 1 Eé 
1 1 
0.0, 6.55 '11:376])-D 
C11 
- Zitg g  +{(0.0261),(-6.55),(-11.876)} 
Dy Can 
ES 
7 2 = Z+ ——— g 
“eq Go Ep» 2 CTi 


EE has three known poles at 0.0261, -6.55, -11.876; the 


other poles were the roots of equation (4-9) 


PT = 0 
SR c 


The zeros of E were the roots of equation (4-8) 


det*G 
[e c D = 0 
Bos C11 


and therefore 


K_(S8+2.3)(S+6.8+j0.8)(S+11.99) 
Gea © (8-0.0261)(3+1.25)(8+2.3)(8+6.55)(S+6.7/5)(S+11.45)(SF11.876) 


and from 
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+ ° 
PA, PNE dC A 


+ 
eq 5 nn Bon: 
0.0394(S+9. 55) 0.144(S-0.02325)(S+6.51)(S+12.02) 
E (500260) (50.55) (S41 S70)NNS(S-0T0261)- CG:6.55) 4. CS:11.876)? 
| 4 .9.08(8*0. 72) (86.52) > 
S(S-0.0261)(S+6.55)(S+11.876) 
P -0.0394(S+Z, )(S+Z, )----- 
OP Seen 
X__ = -0.0394 
eq (4-19) 


G ^ UMOSGO9A4Ot2s99(5t-T6.85]08S)(BGtll.59) 
eq MS UNO Se S F2.25»(5*6.55)(8t60.75)(St11.65) 69: :11,870) 
Equations (4-10) and (4-19) were the equations used in 
designing the compensator x 
Figure 22 is the root locus plot for designing E22 
The positive gain will push the pole at 0.0261 further into 
the right half plane; therefore positive gain cannot be used 
to compensate this system. Negative gain which can make the 
gain of the system equal to 3.881 will set dominant roots at 
-0.7+30.2, which is a pretty good position. Therefore, 
using the gain of the system equal to 3.881, we can make 


system stable and the pitch angle loop has good damping. 


The compensator En thus becomes a pure gain of value 
282 


K2 K = 3.881 
eq 
MS 
K2 = -0.0396 7 798 
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The configuration of the compensated system is shown in 


Figure 23. 





FIGURE 23 


BLOCK DIAGRAM SHOWING THE COMPENSATED SYSTEM 


Referring to Figure 20, the dominant roots were at -1.25 
and -2.2 which gave a, 7 1.658, and z = 1,04, Figure 24 
shows the response of the craft when the height above the 
water has changed from 9.5 ese 17.58feet, and FPigure 25; 
when it changes from 11.5 feet to 9.5 feet. The response of 
the craft was the same in both cases and followed very 
closely to the response of the best second order approxima- 
tone 

Figure 26 is the response of the pitch angle from 2° bow 
up back to 0° reference and Figure 27 is the response from 
2° stern up back to 0° reference. The response in both 
cases was the same but faster than the response of the best 


second order approximation. 
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FIGURE 24 


STEP RESPONSE OF HEIGHT (2 FT. UP) 
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FIGURE 25 


STEP RESPONSE OF HEIGHT (2 FT. DOWN ) 
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STEP RESPONSE OF PITCH ANGLE (BOW UP) 
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FIGURE 27 


STEP RESPONSE OF PITCH ANGLE (STERN UP) 
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C. ROOT LOCUS DESIGN FOR ROLL MODE 


After neglecting the small terms, equation (2-4) becomes 





And from equation (2-41) 


L F L F L L 


L = Foppb vsF "vps-zps "vsstzss Fycslycs 


zpr'vpr'tzsr 


Assume that the roll mode was decoupled from the other 





modes 
L= Frprypr'.zsFlyse 

= [-D2(0.0028A;*0.018)65] Ly 5, * [7D2 (0.002845 *0.018)6 5] Ee 
Sp = 79g 
L - —29959.49gt(-29959. 46.) 

= -59918. 805 
Pp =$ =-0.18, 

Taking the Laplace Transform, the transfer function is 
simply 

$e i -l | (4-20) 
S Š 


Figure 28 is the root locus plot after compensation 
using a lead network which has zero at -1.5 and pole at -195. 
With the gain of -82, the roots were at -3.5:30.75, which 
gave the value of % = 0.96 and De 3.6. Figure 29 is the 


step response roll from 2° port up back to O° reference and 
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Figure 30 from 2° starboard up to 0° reference. The same 
response is obtained in both cases. The block diagram of 


the automatic control of this system is shown in Figure 31. 
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FIGURE 29 


STEP RESPONSE OF ROLL ANGLE (PORT UP) 
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FIGURE 30 


STEP RESPONSE OF ROLL ANGLE (STARBOARD UP) 
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VI. CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER STUDY 


The transfer functions and root locus study has shown 
that this system is unstable by itself because of the open 
loop pole in the right half plane. From studying the effect 
of each foil on the system, it has been shown that stabiliza- 
tion can be achieved by simply closing eacn loop and using 
the appropriate gain. However, that will stabilize only the 
effect of one foil when the other is fixed. 

To stabilize the systems when both forward and aft foils 
are used, the total system has been treated as a multivari- 
able system. To reduce the work in designing the automatic 
control system, some assumptions had to be.made, such as no 
line between the pitch-heave modes, roll mode and yaw 
mode. Thus, the system will have only a 2x2 transfer func- 
tion matrix for the pitch-heave mode, and single loop for 
the roll mode. The step response study has shown that the 
responses of the compensated system followed closely the re- 
sponse of the best second order approximation. The system 
was made stable in both pitch-heave and roll mode. However 
this result cannot insure that the system is completely 
stable in six degrees of freedom. 

Recommendations for further study are: 

1) Study of yaw motion and the craft turning character- 

istic 


2) Wind effects 
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3) 


4) 


5) 


Produce a simulation for the transition period be- 
tween hullborne and foilborne 

Take the dynamic characteristic of the sensors into 
account 

Study the response of the craft in a sinusoidal sea 


and in a random sea. 
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APPENDIX A 


SUMMARY OF COMPLETE EQUATIONS 


EQUATIONS FOR ACCELERATION IN BODY AXIS 


° Ap = 

U = =F_+RV-QW CASIO 

Ÿ = lx +pw-RU (A-2) 
m Y 

W = ZE _+QU-PV (A-3) 
m Z 

p.a n ú : = 

P = ril QRCI,,-I, ,)*(R*QP)I..,] (A-4) 

i $ _rop_f n 


= 


VELOCITY IN BODY AXES 


U = f Udt (A-6) 
V = f Vat (EN 
W = Sf Wat (A-8) 
p = NEP (A-9) 
Q = f Qdt (A-10) 
Ss kat Celis) 


TRANSFORMATION OF VELOCITIES FROM BODY TO EARTH AXES 


Up=Ucos0cosy+V(cosysingsing-sinycosg)+W(cosysinocosg+sinysinò ) 


(A-12) 
V„=Ucosösiny+V(cosycosd+tsinysindsind)+W(sinysindcosd-cosysing) 

(A-13) 
W,=-Usin0+Vcos@sing+WcosAcosy (A-14) 


E 
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P+ÿsine 


De Ge 
il il 


Qcosó-Rsinó 


== 
li 


(Qsinò+Rcosò )cos8+($-P)sin9 


POSITION IN EARTH AXES 


Xp = J U,dt 
vM vot. 
Z, = S W,dt 
d = S dat 

0 - f dat 

y = f pat 


EXPANSION OF FORCE AND MOMENT EQUATIONS 


22 "oe LS x 

e rans ER 

Fa è Fog t "Ez 

E e e ps uz ss Oy 
a CE oz pie das Iza 
N = (Flore 


EXPANSION OF GRAVITY TERMS 
mg, = -mg sino 
mg, = mg cosgsinó 


mg; = mg cosgsing 
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Lz)cs 


(A-15) 
(A-16) 


(A-17) 


CASE 
(A-19) 
(noo 
235) 
(D 


(A-23) 


(A-24) 
(A-25) 
(A-26) 
(A-27) 
(A-28 ) 


(A-29) 


(A-30) 
(A-31) 


(A-32) 





FOIL VELOCITY COMPONENTS IN BODY AXES 


Center Foil: 
Sa cp € 
CA cri cn e 
Me cp S 

Port Foil: 
Apr pri 
Vp = V-La pp PtLypp' È 
PP p: cpi 


Starboard Foil: 


Ug UE gp 8 Lygp'È 
Vs = V-Lygp'P+Lygp'R 
Wo = W+Lygp'P-Lygp*@ 
Mid Foil: , 
Uy = UtLyyp°@ 
Vu = Vokzur **xyp P 
Wu © W-Lyyp°Q 
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(A-33) 
(A-34) 


(A-35) 


(A-36) 
(A-37 ) 


(A-38) 


(A-39) 
(A-40) 


(A-41) 


(A-42) 
(A-43) 


(A-44) 





TRANSFORMATION OF WATER PARTICLE ORBITAL VELOCITY FROM EARTH 


TO BODY AXES 


Uer = Unypcosbcosd+V „„„Ssinbcosd-W sino (A-45) 

Ver = Uny1](cosysin8sinp-sinycosó )+Viy, y (sinysin0sing+cosycosó) 
TWoyrcosesin$ (A-46) 

Rut = U,yy(eosvsindcosd+tsinysind)+V „„„(Sinysindcosd-cosysind) 
+Wypcos9cosd (A-47) 


Each of these equations must be repeated for each foil. 


RELATIVE VELOCITY COMPONENTS 


Uni = u = Ur (A-48) 
VRI = VI = Vor ( A- A9) 
WRI = Wi = WI (A-50) 


One set of equations for each foil. 


ANGLES OF ATTACK AND SIDE SLIP 


= = * 
a arctan WRI/URI * Or FIXED (A-51) 
V 
By = arcsin === (A-52)* 
+ + 
AA 


* One set of equations for each foil. Angles of attack and 


side slip are calculated in radians. 
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TOTAL RELATIVE VELOCITY AT A PARTICULAR FOIL OR STRUT 


V YU "e Wo (A-53) 


HYDRODYNAMIC FORCES IN WATER AXES 
(Each foil equation adn strut equation must be repeated for 


foil and strut, respectively) 


Foil Lift: 
= 1 2 sch 
Fr =p V ACH (A-54) 
Foil Drag: 
= 1pny? — 
For 20V ALC. (A-55) 
Strut Drag: 
“= 1 2 | A 
Fos PV A Cha (A-356) 
Í 
Strut Side Force: 


S S S 


TRANSFORMATION OF FOIL FORCES FROM WATER AXES TO BODY AXES 


Tom = -Fippe0s0cosB+F] jpSina (A-58) 
Bart = -Fy1pSinB (A-59) 
FIF = =F, ¡pSinacosg=F, ¡peosa (A-60) 
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APPENDIX B 


PLOTS OF HYDRODYNAMIC COEFFICIENTS 


Fwd. Foil Lift Coefficient vs. Angle of Attack 
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Mid Foil Lift Coefficient vs. Angle of Attack 
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Aft. Foil Lift Coefficient vs. Angle of Attack 
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Forward Foil Drag Coefficient vs. Angle of Attack 
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